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Abstract 

We demonstrate that the valence energy-loss function of hexagonal boron nitride (/iBN) displays 
^ • a strong anisotropy in shape, excitation energy and dispersion for momentum transfer q parallel 

or perpendicular to the /iBN layers. This is manifested by e.g. an energy shift of 0.7 eV that 
cannot be captured by single-particle approaches and is a demonstration of a strong anisotropy 
in the two-body electron-hole interaction. Furthermore, for in-plane directions of q we observe a 
splitting of the vr-plasmon in the FM direction that is absent in the FK direction and this can be 
traced back to band-structure effects. 



Layered hexagonal boron nitride (/iBN) is the III-V compound counterpart of graphite. 
The recent interest in hBN is partly due to Watanabe et al. [1] who showed that hBN exhibits 
the potential for lasing at high energies (5.8 eV). This makes hBN an attractive candidate 
for optoelectronic applications in the ultraviolet energy range. Furthermore, a single layer 
of /iBN is isoelectronic to graphene and can be considered the precursor material for boron 
nitride nanotubes. Despite the relatively simple crystal structure, hBN appears to be a 
challenge to both experiments and theory. The band-structure of hBN displays an indirect 
gap J2|, yj and the optical absorption spectrum is dominated by correlation effects leading 
to a strong Frenkel-type excitonic peak at 5.8 eV |3H5| in agreement with the experimental 

le absorption and luminescence spectra 



-|8|]. Recent luminescence experiments 



findings (l|. The origin of the fine structure of t 

around this peak, however, is still under debate 

point towa.ds the .ole of defect. M - agree„.e„t w.th the theoretical .ugge.tion .n 

ref . 16|. 

A major source of experimental difficulties lies in the growth of large defect-free hBN 
crystals. This is especially important for the electronic properties since hBN stacking order 



can have a substantial effect on the band structure |12l . Il3| . In the case of cubic boron 
nitride, the discrepancies between experimental and theoretical results can often be traced 
back to the use of low quality samples [ij] . 

In order to study the dynamics (e.g. possible dispersion) of the valence excitations one 
needs to use inelastic scattering probes such as electron energy loss spectroscopy (EELS) 
or non-resonant inelastic x-ray scattering (NRIXS). NRIXS has the advantage to study the 
dispersion of low lying excitations beyond the first Brillouin zone (1st BZ). This has been 
used numerous times in the past (see e.g. [l5|, [l6|) and gives new important information on 



the properties of excitation, at sub-u,nt-cell length .c.le^. When studying ,ong.tnd.na, 

excitations one often finds non-parabolic dispersions [18[ and periodicity for low energy 



plasmons 19|, l20[. The valence bands and the lowest conduction bands of hBN can be 
understood as combinations of the a and vr states of the hexagonal boron nitride sheet. 
At low momentum transfer the low-energy-transfer structures of the loss function can be 
divided into transitions between states with same parity {a-cx*, tt-tt*) when the momentum 
transfer is in the plane and different parity (vr-cr*, a-ir*) for momentum transfer along 
the c-axis 2l|. The strong anisotropy in the electronic response can also be seen in the 
difference for the dielectric constants ^oo parallel and perpendicular to the planes (4.40 and 



2.53, respectively |3[). 

In this work, we study the loss function of /iBN using NRIXS. We show detailed experi- 
mental results of the dispersion of various features and plasmons for the momentum transfer 
q in different crystallographic directions within and beyond the 1st BZ. The loss features 
show a strong directional dependence not only in the comparison in-plane/out-of-plane but 
also for the different directions within the plane. The origin of the different spectral features 
and their direction and momentum dependence are analyzed by ab-initio calculations at the 
level of the random-phase approximation (RPA) and with methods of many-body perturba- 
tion theory. In order to exclude a perturbation of our results by the low crystal quality, we 
have characterized the sample in detail. 

The sample was a colorless and transparent hexagonal single crystal /iBN platelet about 
700 yum wide and 70 /im thick. The FA direction was found to be very nearly parallel 
to the normal of the platelet. We used single-crystal X-ray diffraction to verify that the 
sample exhibited an AB-type stacking, thus ruling out other stable or metastable stacking 



sequences 
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13| 



The NRIXS spectra were measured on ID16 at the European Synchrotron Radiation 
Facility, Grenoble, France. The experiment was carried out using the eV-resolution spec- 
trometer 



22 



23l |. An energy resolution of 0.6 eV was determined from the FWHM of the 



quasi-elastic line. The measurements were performed using the inverse scan technique 24 1. 
The momentum transfer resolution was Aq/q ^0.17 near the K point. 

Theoretical NRIXS spectra have been calculated at the level of the random-phase ap- 
proximation (RPA) starting with the Kohn-Sham DFT wave-functions in the local-density 
approximation (LDA) [25,]. The dielectric function is obtained through eG,G'(q;i^) = 
1 — f (q + G)xg g'(^' '-^)' where f (q) = 47r/|qp is the Coulomb potential in reciprocal space 
and Xg g'(^5 ^) i^ ^^^ independent particle polarizability which is obtained from a sum over 
transitions from occupied to unoccupied bands [26|. G, G' denote reciprocal lattice vectors, 
and q is restricted to the 1st BZ. 

The loss function (T(q, w) is then obtained from the inverse of the dielectric matrix 



a(q+ G,a;) = -lm{e(.]G{^, uj)) , 



(1) 



where q = q+G. Crystal local- field effects are automatica 
the off-diagonal elements of eG,G' in the matrix inversion 



ly included by taking into account 
27| . In order to compare with the 
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FIG. 1: (Color online) The experimental (black dots) and theoretical (red lines) NRIXS spectra in 
the 1st BZ along three crystallographic directions. The g^-weighed spectra are vertically displaced 
proportionally to the value of momentum transfer, which is indicated on the vertical axis. The 
dashed lines are guides for the eye for features discussed in the text. The inset shows the 1st BZ 
of /iBN with the irreducible part shaded. 



experimental loss function, we have used a broadening rj = OA eV. In silicon the use of the 
time-dependent LDA kernel leads to a considerable enhancement in the description of short- 



range exchange- correlation effects 28[ and to a better agreement with the experiments. We 
have checked that the use of time-dependent LDA does not give a substantial improvement 
in the case of /iBN . 

The experimental and theoretical NRIXS spectra along the three crystallographic direc- 
tions within the first are depicted in Fig. [H The experimental spectra were normalized 
to have the same area as the calculated spectra in the energy range shown in the figure. 
Overall the calculations reproduce the experimental spectra well. The ab-initio calculations 
match the experimental feature positions, their relative weights as well as the momentum 
transfer dependence. From the figure it is evident that the spectral features along the FA 
direction do not exhibit any significant dispersion. The high anisotropy of /iBN is clearly 
reflected in the differences between the spectra in the hexagonal plane (FK, FM) and the 
spectrum perpendicular to it (FA). For q ^ 0.6A~^, the spectra along FK and FM are nearly 
identical. However, as the value of q is increased an anisotropy also within the hexagonal 
plane is clearly observed. In the high-energy range this anisotropy shows up mainly as a 



different rate of relative spectral weight increase for the feature between 35 eV-40 eV. The 
in-plane anisotropy is most evident in the behavior of the vr plasmon. Its energy disperses 
from 9 eV at small values of momentum transfer to about 12 eV when q is near the boundary 
of the 1st BZ. Along FM an additional peak around 8 eV develops for q > 0.87 A^^ while 
in the FK direction only a weak shoulder is detected. The double-peak structure along FM 
is also visible in Im(e(q)). This indicates that the directional anisotropy can be interpreted 
in terms of interband transitions. Fig. [2] (a) shows the band structure of hBN. For q = TK, 
the transitions that dominate the plasmon at 12 eV are the ones from the vr band at A/H to 
the vr* band at H/A, respectively (red arrows). For q = FM, the dominant transitions are 
the ones from the vr band at A/L to the vr* band at L/A (blue arrows). The observed 8 eV 
transition at q = FM originates from a vr — tt* interband transition from L to a neighbor- 
ing high symmetry point L'. This is marked by the vertical green arrow in Fig. [2] (a) even 
though it is not a vertical transition but one with with a momentum difference of q = FM. 
At q = TK the plasmon peak around 8 eV is missing because the joint-density of states 
displays only a minor peak there. 

The FM and FK theoretical spectra have been blueshifted by 0.8 eV whereas for the 
FA spectra a blueshift of 1.5 eV was used. These shifts reflect the well-known fact that 
the DFT band-structure usually underestimates the transition energies between occupied 
and unoccupied bands. For /iBN it was shown j^, y, |6| that electron-electron correlation 
(calculated on the level of the GW-approximation) increases the transition energies by about 
2 eV with respect to DFT-LDA calculations. At the same time, the attractive electron-hole 
interaction reduces the transition energy such that the difference between experimental and 
theoretical NRIXS spectra is less than 2 eV. We note that the GW-correction alone cannot 
explain the observed anisotropic shift since vr and a bands are renormalized in the same way. 
Therefore the experimental NRIXS data cannot be explained with single particle theories 
and the explicit inclusion of electron-hole interaction is necessary. 

In order to check if the combined effect of electron-electron and electron-hole interaction 
does indeed explain the (anisotropic) shift of the RPA-spectra, we carried out calculations 
at selected momentum transfers including excitonic effects on the level of the Bethe-Salpeter 



Equation (BSE) [29|]. We used the approach of Refs. 



|30| and converged a 
parameters [3l!]- To approximate the quasi-particle band energies of Ref. 
"scissor" of 1.92 eV to shift the LDA conduction bands energies and a small stretch of 5% 
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FIG. 2: (Color online) a) DFT-LDA band-structure of /iBN. On the right the 2D projections of 
the electron probability density |^^(r/i,re)p for the lowest bright exciton in the optical absorption 
spectrum are shown with q ^ b) parallel and c) perpendicular to the /iBN planes. The hole 
position r/i (marked by white circle and arrow) is fixed 0.4 a.u. above a nitrogen atom. Balls and 
sticks indicate the atomic layers. Calculations of the excitonic wave functions have been performed 
using BSE. 

was applied to the valence bands. As shown in Fig. |3] the BSE results without additional 
energy shifts agree well with the RPA spectra that were blueshifted to the experimental 
energy scale. The slight differences between the two calculations seem to originate mainly 
from different weights of various spectral features. Also, it should be noted, that (except for 
the energy shift) no extra features seem to arise in the BSE results in comparison with the 
RPA calculations. This confirms that the RPA successfully describes long-range excitations 
even in strongly anisotropic and layered materials like /iBN, with the BSE leading only 
to minor corrections. Electron-hole effects can alter extended excitations only when the 
attraction is very strong, like in a wide-gap insulators such as LiF 29 1. 

We found that the layered structure of hBN leads to a strongly anisotropic electron-hole 
interaction. This is demonstrated in Fig. [2] b) and c) where we show the excitonic "wave 
functions" for the lowest visible excitations in the optical absorption spectrum. For q parallel 
to the planes (gy), electron and hole are localized within one plane forming a strongly bound 
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FIG. 3: (Color online) The comparison of the current RPA calculations (black solid curve) with 
a BSE calculation (red dashed curve) at selected momentum transfers. The RPA spectra are 
blueshifted as discussed in the text. Experimental data is also shown (black dots). 

quasi-two dimensional compound |3|, |6[. For q perpendicular (gx), electron and hole are 
localized on different layers and form a more weakly bound three-dimensional compound. 
This leads to a larger exciton binding energy for gy than for g_|_ and thus explains why the 
upshift due to the combined e-e and e-h interaction is lower for q\\ than for q^. 

The NRIXS spectra with a larger range of momentum transfers directed along FK are 
displayed in Fig. HI For most parts the calculated spectra agree well with the experimental 
ones. The prominent peak dispersing between 25-30 eV in the experimental spectra is 
reproduced very well by the calculations throughout the measured range of momentum 
transfers. At high momentum transfer (> 2 A~^) a new peak showing a small dispersion 
appears just above 20 eV. Its appearance and dispersion is reproduced well by the calculation. 
In the same momentum transfer range calculation shows a peak at around 15 eV which in the 
experiments appears to be less well pronounced and at a few eVs higher in energy. Around 
7 eV in the high momentum transfer experimental spectra a sharp peak is observed, which 
is not present in the calculated spectra. 

In conclusion, we have studied the anisotropy in the valence electron dynamics of /iBN 
using NRIXS. Besides the in-plane/out-of-plane anisotropy, the in-plane spectra show a 
directional dependence at large q. In particular, the vr-plasmon is split along the high 
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FIG. 4: (Color online) The theoretical (left) and experimental (right) NRIXS spectra with the 
momentum transfer directed along the FK direction using several values of q. 

symmetry line F-M but not along the line F-K. The anisotropic splitting is an effect of the 
anisotropic dispersion of the vr-bands. We have detected an anisotropy in the electron-hole 
interaction for momentum transfer parallel/perpendicular to the layers. For gy, the electron 
and the hole tend to form a strongly bound quasi-two dimensional compound while for q^ 
electron and hole tend to be localized on different layers and form a weakly bound three- 
dimensional compound. This anisotropy explains the 0.7 eV difference in the blueshift that 
must be added to the RPA spectra (which neglect the effect of electron-correlation) in order 
to match the experimental data. 
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